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Abstract 

Spectral statistics in band structures is studied in a realistic model describing su- 
perconducting quantum interference devices (SQUID's). By controlling an external 
magnetic flux, the level statistics may show a crossover from the GUE to the GOE. 
Effects of secondary discrete symmetries seen in specific regions of the first Brillouin 
zone are also discussed. 
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Recent development in physics of quantum chaos has its origin in studies of 
spectral statistics of quantum systems whose corresponding classical dynam- 
ics exhibits chaos. It is now well known that statistical properties of energy 
spectra are universal in various quantum systems including, e.g., atoms, nu- 
clei, and electron billiards as well as classical wave-mechanical systems such as 
microwave cavities [1]. Recently, spectral statistics in band structures in ex- 
tended periodic systems has been studied for driven systems [2] , electronic sys- 
tems [3,4,5,6] and photonic crystals [7]. Band structures have also been studied 
in the context of avoided band crossings [8]. Most of these band-structured 
systems, however, are sensitive to the accompanying impurities and/or inho- 
mogeneity. It is highly desirable to find corresponding systems robust against 
the extrinsic randomness. 

In this Letter, we consider superconducting quantum interference devices 
(SQUID's), whose periodic features are described by a few macroscopic vari- 
ables called superconducting phases. We here focus on small SQUID's with 
three junctions where a two-dimensional periodic potential for the phases can 
be realized. Because of the macroscopic nature of the phases, the potential is 
robust against impurities and/or inhomogeneities of the system. This feature 
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Fig. 1. (a) The three-junction system considered in this Letter. The states of the 
SQUID are controlled by a magnetic flux $ext and gate voltages, V^'s. The flux 
is vertical to the page, (b) The equivalent circuit of the SQUID. The Josephson 
energy and capacitance of the i-th junction (shown by x) are denoted by Eji and 
Ci, respectively. The two upper superconducting islands are coupled to the voltage 
gates through capacitances, Cgi and Cg2, respectively. 

is of great advantage to study the spectral statistics; This is particularly in 
contrast to the two-dimensional electron systems where impurity scattering 
and inhomogeneity are unavoidable. The SQUID system also has the advan- 
tage of controllability. Potential shapes for the phases can be controlled by an 
external magnetic flux, and the Bloch wave-numbers can be chosen by voltage 
gates. In the limited context of application to quantum computations, one 
is concerned only with the lowest few quantum levels [10,11,12,13]. Moreover 
generally, however, if level spectroscopy at a high-energy region will become 
possible, this system would give a good stage for studies of quantum chaos in 
band structures. 



Already we investigated the classical phase dynamics of SQUID's with two 
junctions and discovered the mysterious normal diffusion phenomena there: 
The diffusion coefficient was found to take a fractal-like feature [14]. SQUID's 
with two junctions, however, were on scale of lOOyum and thereby classical. 
On the contrary, SQUID's with three junctions considered in this Letter is on 
scale of 1/im and their quantization is highly desirable to elucidate a quantum 
signature of the underlying mysterious normal diffusion phenomena. 

The configuration of the SQUID system is shown in Fig. 1 (a), while the equiv- 
alent circuit is shown in Fig. 1 (b). The loop of superconductor includes three 
Josephson junctions denoted with x in the circuit. The Josephson energy and 
capacitance of the i-th junction are denoted with Eji and Ci, respectively. The 
two upper superconducting islands are coupled to the voltage gates through 
the capacitances, Cgi and Cg2, respectively. The states of this SQUID are con- 
trolled by an external flux penetrating through the loop, ^ext and by gate 
voltages, V^'s. 

By controlling fabrication conditions, both the Josephson energies and capac- 
itances can be varied. In this Letter, for simplicity we assume Eji = Ej2 = 
Ej3{= Ej). We also assume 6*1,6*2,6*3 <C C*gi = C*g2(= C*g), where the junc- 
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tion capacitances can be neglected. These assumptions are not essential to the 
discussion of the classical-quantum correspondence through spectral statistics. 
We assume these conditions because it is convenient to study effects of discrete 
symmetries such as parity. 

Behaviors of the Josephson junction are described by the superconducting 
phase difference at the junction. We denote the phase difference at each junc- 
tion with 0j. The directions to define these differences are shown by the arrows 
in Fig. 1 (a). Then, the Hamiltonian of this SQUID consists of three parts: 



where Hq, Hj and i^L are a charging energy, a junction energy and an induc- 
tance energy, respectively. First, we show that the latter two parts correspond 
to a potential energy for the phase differences. The junction energy is deter- 
mined by the phase differences as 



while the inductance energy is determined by the persistent current along the 
loop, and given as 



Here, $o = ^/2e is a unit flux, and L is an inductance determined by the size 
and shape of the loop. Here, we consider a small loop with a small inductance 
satisfying ^l/2L ^ Ej. Such SQUID can be easily realized experimentally by 
fabricating the superconducting loop with the size of order of micrometers. 
In this small SQUID, the inductance energy approximately gives a constraint 
01 ~ 02 ~ 03 ~ 27r$cxt/'^'o = on the phases. By using this constraint, the 
potential energy for the phases are obtained as 

Hj = -Ej cos 01 - Ej cos 02 - Ej cos(0i - 02 - 27r/) , (4) 

where / = $ext/^o is a normalized external flux. This two-dimensional peri- 
odic potential has a square unit cell defined by < 0i < 27r and < 02 < 27r. 

On the other hand, the charging energy He is described by the charges of 
the junctions. We denote the charge induced at the i-th gate capacitance by 
Qi = —2eni, where rii is the number of the induced Cooper-pair. It is known 
that the number of the Cooper pairs, rii is conjugate to the phase difference 



H — He + Hj + Hi^, 



(1) 



Hj = -Ej cos 01 - Ej cos 02 - Ej cos 03, 



(2) 




(3) 
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0i and the following exchange relation satisfies: 

[(pi,nj] = iSij, (l<i,i<2). 



(5) 



Since the junction capacitances are neglected, the charging energy is calculated 
as [10] 

He = 2^(gi + Ql) + V,Q, + V2Q2. (6) 

This energy is rewritten by using the numbers of Cooper-pairs as 

He = 4Ec{n, - nlf + 4Ec{n2 - nlf, (7) 

where Eq = e^/2Cg, and n* = CgVi/2e is a normalized gate voltage. This 
gives a kinetic energy for the phases. 

Here, note that the gate voltages can control the Bloch wave-numbers. This 
can be checked by considering the wave function ^2) satisfying the peri- 

odic condition for both (f>i and (f>2. For the alternative wave function ^2) 
defined by ^2) = e-^"^'^i-^"2'^2^(0i, (j)2), the kinetic part of the Hamilto- 

nian for ^2) is modified as 

He ^ ^Ecnl + 4Ecnl (8) 

while the wave function .^(0i,02) satisfies the following Bloch theorem 

C{<Pi + 27rmi, 02 + 27rm2) = e'^'^'^+'"''^H{<Pi, h)- (9) 

Hence, the normalized voltages, n*'s correspond to the Bloch wave-numbers 
for the Hamiltonian H = Hq + Hj. We find that the feature of classical phase 
dynamics for H leads to the same mysterious normal diffusion phenomena 
in squid's with two junctions [14]. (We suppress the classical result in this 
Letter.) 

The energy spectrum of the SQUID is obtained by diagonalizing the Hamil- 
tonian H — He + Hj for fixed values of n* and 712- We adopt the charge 
eigenstates as the base wavef unctions. Then, it is convenient to rewrite the 
Hamiltonian by introducing an annihilation operator bi — e^'^' satisfying an 
exchange relation [rii, bj] = —bidij as 



H^4Ec{ni-niy + 4Ec{n2-niy 

- f (&i + b{) - f (&2 + bl) - ^{bAe-'^'' + b2b\e'^-^f). (10) 
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The charging energy gives diagonal parts, while the junction energy off-diagonal 
parts. Here, we should note that if / is an integer or half-integer, the Hamil- 
tonian can be represented by a real matrix. 

In this Letter, we take Eq = 0.002£'j, and consider the case of / = 0.5 
and / = 0.25. The classical phase dynamics based on Eq. (1) shows chaotic 
behaviors in the dominant region of the corresponding Poincare section, when 
the energy is taken between the saddle-point potential energy i^sad and the 
potential-top energy E'top- Here, the spectral properties are studied in the range 
Egad ^ E < Egad + 0.75(£'top — -E'sad)- The ordinary unfolding procedure is 
performed to obtain the level statistics [15,16]. Depending on the Bloch wave- 
numbers, nl and the energy spectra show various statistics: the Poisson- 
like, GOE, and GUE statistics. To see this wave-number dependence, we fit 
the spectra to the Izrailev distribution [17] 



PiS) = A{—j exp 



which interpolates among different universality classes by the parameter (3. 
The parameters, A and B are normalization constants. One can recover the 
Poisson distribution for P — 0, while it agrees with the GOE (GUE) distribu- 
tion for f3 — 1 — 2) on the 5 percent level. Although the Izrailev distribution 
is originally proposed for localization problems, we here use it for qualitative 
characterization of the spectra. 

The density plots of the fitted values of /3 in the first Brillouin zone (FEZ) of 
the 71-1-^2 plane are shown in Fig. 2. In the case of / = 0.25, the spectra shows 
the GUE distribution {(3 = 2) except for the symmetry lines, n* = ^2 and nl = 
-n*2, and the points, (n^, nl) = (±0.5, 0) and (0, ±0.5). In the case of / = 0.5, 
the spectra shows the GOE distribution {(3 = 1) except the same symmetry 
lines and points. Let us first discuss the universality class of the region, which 
is on neither the symmetry lines nor points. In the case of / = 0.25, the 
Hamiltonian includes the imaginary number in the last term in (10). Hence, it 
is expected that the time-reversal symmetry is broken, and the level statistics 
is governed by the GUE. Actually, the nearest neighbor distribution P{S) and 
spectral rigidity A3(L) agree well with the GUE statistics as shown in Fig. 3. 
In the case / = 0.5, the Hamiltonian can be expressed by a real symmetric 
matrix as seen in (10). Hence, the system exhibits the time-reversal symmetry, 
and the level statistics belongs to the GOE. Actually, the result for P{S) and 
A3(L) agrees with the GOE statistics as shown in Fig. 4. Since the value of 
/ can be varied by the external magnetic fiux, the universality class can be 
easily controlled in this SQUID system. 

Thus, the signature of quantum chaos has been obtained in the SQUID system. 
So, the main purpose of this Letter has been achieved. It, however, remains 
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Fig. 2. The density plot of the parameter /? obtained by fitting the spectrum to the 
Izrailev distribution: (a) / = 0.25 and (b) / = 0.5. /? = (the darkest) and (5 = 2 
(the brightest). 

to explain the deviations from the expected universal class on the symmetry 
lines and points as seen in Fig. 2. We end this Letter by brief discussion on 
this secondary feature in the SQUID system. 

First, we consider the case of / = 0.25. We divide the regions where the de- 
viation from the GUE statistics is seen into five parts: (A) the line n\ = 
(n* ^ 0), (B) the line nl = -n^ (n* ^ 0), (C) the point {n\,n1^ = (0,0), (D) 
the points (n^,^;) = (±0.5,0), and (E) the points (nj,^) = (0,±0.5). The 
deviation from the GUE statistics in each region can be explained by corre- 
sponding additional symmetries characterized by operators commuting with 
the Hamiltonian. In the region (A), it is crucial to consider the operator Pi 
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Fig. 3. (a) The nearest-neighbor level distribution P{S) and (b) the spectral rigidity 
A3(L) for / = 0.25. The normalized gate voltages are taken as (?^^,?^2) = (0.25,0). 

defined by \n2,ni) = Pi\ni,n2), where \ni,n2) is the eigenkets of the Cooper- 
pair-number operator. It can be checked that the operator PiT commutes with 
the Hamiltonian on the fine nl = where T is the time-reversal operator. 
Then, the spectral statistics is governed by the GOE statistics by the discus- 
sion of the false time-reversal symmetry [18]. Similarly, in the region (B), the 
mirror- reversal operator P2 defined by | — ^2, — ^^i) = -P2|^i; ''^2) commutes with 
the Hamiltonian. Then, the eigenkets can be classified into two kinds by the 
parity defined as the eigenvalue of the operator P2. The spectrum taken only 
for even(odd) states is expected to show the level repulsion. The level correla- 
tion between odd and even states is, however, suppressed because the matrix 
elements of the Hamiltonian {ipilHlipj) is zero when -0, and ipj have different 
parity. As a result, the spectrum become the mixture of two independent GUE 
statistics. This mixed level statistics is denoted in this Letter with 2-GUE [5,6]. 
In the region (C), both PiT and P2 commutes with the Hamiltonian. Then, 
the spectrum becomes the mixture of the two independent GOE statistics 
(2-GOE). In the region (D), a kind of the mirror reversal symmetries exists. 
The symmetry operator P3 is defined as |1 — ni, —71-2) = P3\ni,n2), and P3T 
commutes with the Hamiltonian. Then, the spectral statistics obeys the GOE 
by the same reason as the region (A). Similarly, in the region (E), P4T com- 
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Fig. 4. (a) The nearest-neighbor level distribution P{S) and (b) the spectral rigidity 
A3(L) for / = 0.5. The normalized gate voltages are taken as (n|,n2) = (0.25,0). 

mutes with the Hamiltonian, where P4 is defined as | — ni, 1 — = -Pil^i, ^2). 
Hence, the level statistics is described by the GOE. 

Next, we consider the case of / = 0.5. The deviation from the GOE statistics 
is observed near the same lines and points as the case of / = 0.25. The 
explanation of this deviation is quite parallel to the previous paragraph. In 
the region (A) and (B), the operator Pi and P2 commute with the Hamiltonian, 
respectively. Hence in these regions, the level spectrum becomes the mixture 
of two independent GOE statistics (2-GOE). In the region (C), because both 
Pi and P2 commute with the Hamiltonian, the spectrum obeys the 4-GOE 
statistics. In the region (D) and (E), the spectrum is described by the 2-GOE 
statistics, since P3 and P4 commute with the Hamiltonian, respectively. 

In summary, spectral statistics in band structures has been discussed in a 
realistic SQUID model. This system can change the Bloch wave-numbers by 
controlling the gate voltages. The normalized external flux / can control the 
time-reversal symmetry of the system, and the crossover from the GUE to 
the GOE can be observed by varying /. By taking the circuit parameters 
as symmetric, the effect of the secondary discrete symmetries can also be dis- 
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cussed. On specific lines and points in the first Brillouin zone, the spectra may 
obey different statistics determined by the operator characterizing the discrete 
symmetry of the system. We expect that this work can be a starting point to 
consider the quantum chaos in a well-controlled system. A variety of dynami- 
cal features of quantum chaos can also be observed in the quantum dynamics 
by using the excellent ability of quantum manipulation and measurement in 
the SQUID system, which will be investigated in due course. 

The authors thank A. Terai for useful discussion and critical reading of the 
manuscript. 
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